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Abstract 

We study the Voronoi and void statistics of super-homogeneous (or hyper uniform) point patterns 
in which the infinite-wavelength density fluctuations vanish. Super-homogeneous or hyperuniform 
point patterns arise in one-component plasmas, primordial density fluctuations in the Universe, 
and in jammed hard-particle packings. We specifically analyze a certain one-dimensional model by 
studying size fluctuations and correlations of the associated Voronoi cells. We derive exact results 
for the complete joint statistics of the size of two Voronoi cells. We also provide a sum rule that 
the correlation matrix for the Voronoi cells must obey in any space dimension. In contrast to 
the conventional picture of super-homogeneous systems, we show that infinitely large Voronoi cells 
or voids can exist in super-homogeneous point processes in any dimension. We also present two 
heuristic conditions to identify and classify any super-homogeneous point process in terms of the 
asymptotic behavior of the void size distribution. 



I. INTRODUCTION 



Point patterns are ubiquitous in Nature. Examples include those defined by the coordi- 
nates of the particles in a many-particle system, such as the molecules of a liquid or crystal, 
stars of a galaxy, or trees in a forest. Understanding how the number of points fluctuates at 
a given length scale reveals important structural information about the point pattern. Such 
local density fluctuations have been studied for a variety of physical systems, including one- 
component plasmas molecular liquids |2], and the large-scale structure of the universe 

Point patterns in which the infinite-wavelength density fluctuations vanish, referred to 
as super-homogeneous B or hyperuniform Q], are of particular interest to us in the present 
paper. Regular lattices of points in space are the simplest examples of super-homogeneous 
point patterns, but such point processes are neither statistically spatially stationary (ho- 
mogeneous) nor isotropic. Stochastic super-homogeneous point processes and fluctuations 
have been demonstrated to be verv important in a variety of physical contexts, including the 
study of one component plasmas M , the evolution of primordial matter density fluctuations 
in cosmology y, and the structure properties of jammed configurations of hard spheres sys- 
tems 0]. It is considerably more difficult to construct point patterns that are statistically 
stationary and isotropic, although some examples have been identified . In order to 

shed light on this problem, our general goal is to understand the statistics of the underly- 
ing Voronoi cells associated with the points of stationary and isotropic super-homogeneous 
point processes in arbitrary space dimension d. A Voronoi cell associated with a given point 
consists of the region of space closer to this point than to any other point 

A first step toward the stated goal is to start by examining stationary super-homogeneous 
point processes in one dimension (where isotropy is not an issue). Specifically, we analyze 
such a particular one-dimensional model by studying size fluctuations and correlations of 
the associated Voronoi cells. We derive exact results for the complete joint statistics of the 
size of two Voronoi cells. It is additionally shown that infinitely large Voronoi cells can exist 
in super-homogeneous point processes in one dimension. We also provide a sum rule that 
the correlation matrix for the Voronoi cells must obey in any space dimension. 
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II. PRELIMINARIES 



Before discussing the details of the model, we recall some general definitions of basic 
quantities that are used to statistically characterize point processes (for rigorous definitions 
and analysis see 

A single realization of a point process is completely determined by the stochastic micro- 
scopic density function n{x), which in d dimension, can be expressed as 

n(x) = ^5(x-x,), (1) 

i 

where (5(x) is the usual c?— dimensional Dirac delta function, Xi is the position of the i*^ 
point in the system and the sum is over all of the points. The microscopic density has the 
following integral property: 

I d'^xn{^) = N[V] , 

where V is any measurable set of the space (i.e., the one- dimensional line in the one- 
dimensional case of interest) and A^[V] is the number of points (particles centers) contained 
in that set. 

The statistics of a point process is completely determined by the infinite set of correlation 
functions: 

/m(xi, ...,Xm) = (ra(xi)...n(x„)) , 

for any integer m > 1, and where (...) indicates the ensemble average over all the possible 
realizations of the point process. For definitions of more general correlation functions see 
Ref. |6|]. Clearly, all of the functions are invariant under any permutations of the variables 
Xi, ...,Xm. For n = 1, we have that 

/i(x) = n(x) = (n(x)) 

gives the local average density of points at the spatial position x and characterize completely 
all the one-point statistical properties of the system. However, very often a constant global 
average density is also evaluated through a volume average: 

77-0= lim — [ d'^xn(x), (2) 

which gives the average density of particles in the system as a whole, and where V is for 
example a spherical volume.. Systems in which volume averages as in Eq. ^ are equal to 
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the relative ensemble averages are called ergodic systems. The quantity /2(x, y)d'^xdf'y gives 
the joint a priori probability of finding a point in the volume element d'^x around x and at 
the same time another in the element d'^y around y. It is the most commonly used function 
to study the correlation properties of an empirical particle distribution. 

If all the /„( Xi, ...,Xm) are invariant under a constant translation of all the points, i.e. if 

X-m ) Im 

(xi + xo, ...,Xm + xo) for any xq and m, the stochastic point process is 
said to be spatially statistically stationary (or statistically homogeneous). In most of what 
follows we will limit our considerations to this class of point process. In such instances, 
n(x) = rio > (the condition > excludes fractal point distributions) does not depend on 
X, and /2(x, y) = /2(x — y) depends only on the displacement vector. If moreover the system 
is statistically isotropic I2 depends only on the scalar distance |x — y|. 

A (i— dimensional point-process is said to be ergodic if, for any function 
-F[n(xi), n(x2), n(xi)] of the microscopic density n(x) in the arbitrary points Xi, X2, x/ 
(where / is finite but arbitrary), the following relation holds: 



It is clear from Eq. Q that spatial stationarity is a necessary condition for ergodicity [6|,|8|,|9|. 
Ergodicity is often supposed apriori as a valid working hypothesis in the analysis of (spatially 
or temporally) stationary stochastic processes. 

In order to measure the density-fiuct nation correlations 5n(x) = n(x) — hq between two 
different points in a statistically stationary point process, the covariance function (also called 
reduced two-point correlation function) C (x) is introduced via 



It is simple to show, from that C(x) can be written as 

C(x) = no5(x) + 'n,g/i(x) , 

where no5(x) is the diagonal part of C(x) present in any stochastic point process inde- 
pendently of the correlations between different spatial points and due only to the discrete 
nature of the massive point-particle distribution, while no/i(x), meaningful for x 7^ 0, is the 
non- diagonal part characterizing the real correlation between different points and vanishing 
for X = Ixl ^ 00. The function /i(x) is referred to as the total correlation function in the 

n 

theory of liquids |10|] . 



lim —[ d'^Xo-F[ra(xo + Xi),n(xo + X2), ...,n(xo + Xi)] = (F[n(xi),ra(x2), ...,n(x;)]) . (3) 





(4) 
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Another important quantity, characterizing the relative weight of each Fourier mode to a 
reahzation of the stochastic point process, is the so-called power spectrum s(k) (proportional 



to the so-called structure factor and called also Bartlett spectrum [l^), which is defined 

by 

.(k)= hrn (|5„(k;L)p) , (5) 

where 

-ikx 



is the Fourier element of the density contrast 5n(x) in a cubic volume of size L. It is simple 
to show that if the point process is spatially stationary then s(k) is simply the Fourier 
transform of C(x): 

s(k) = no + nl J d^x h{x.) e''^'"^ = uq + nlh(k) , 



where h(k) is the Fourier transform in t 



le infinite volume of h{x). This result implies the 



so-called Wiener-Khinchtine theorem 0, which states that the covariance function of a 
stationary point process has a positive Fourier transform converging to Hq for sufficiently 
large k and integrable around k = 0. 

Finally if the system is also statistically isotropic also s(k) depends only on k = |k|. 

III. SUPER-HOMOGENEOUS (HYPERUNIFORM) POINT PROCESSES 

Here we briefly review definitions and basic properties of super-homogeneous (or hype- 
runiform) point processes. Given a spatially stationary point process in d dimensions, we 
can define the variance in the number of points in a sphere Q{R) of radius R (the origin of 
the sphere is arbitrary because of the spatial stationarity) as 



a\R) = {N\R)) - {N{R)y , (6) 

where 

N{R) = [ rf'^xn(x), 
Jn{R) 

is the number of point in the sphere Q{R), which is a stochastic function. 

It is simple to show that Eq. (jHI) can be expressed in terms of the covariance function 
C(x) as follows: 

a\R)= f I d'^xd'^yCU-y) (7) 

Jn{R) JQ.(R) 
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Equivalently, we can express the same quantity in terms of the power spectrum s(k) 0| in 
the following way 

a\R) = -^Jd''k\w{k;R)\'s{k), (8) 
where the integral is over all the k space, and 

w{k;R) = f d'^xe-'^-'' 
JniR) 

is the so-called window function. 

All stationary pointprocesses can be classified in terms of the scaling behavior of 
for large R as follows j3[: 

• If 

J d'^xC(x) = s(0) = A>0, 
i.e., the correlations are mainly positive and short ranged, then 

a\R) ~ R"" 

for sufficiently large R (i.e., for R larger than the range of correlations). The proto- 



typical example of this class of systems is the so-called Poisson point process 
which can be generated by randomly placing points in the space with a given average 
density no > in an uncorrelated manner. In this case, it is simple to show that sim- 
ply C(x) = no5(x) and s(k) = riQ. For this reason we call this class of point patterns 
essentially Poissonian. This is the most common behavior for the number fluctuations 
for homogeneous systems in thermal equilibrium (e.g., an ordinary gas in equilibrium 
at high temperature or a liquid away from critical points). 

If, instead, 

j d'^xC{^) = s(0) = +00, 
with s(k) ~ k~'^ for sufficiently small k where, for deflniteness, < 7 < ci, then 

a\R) ~ R''+^ 

for sufficiently large R. In this case, two-point correlations are again mainly positive 
but are long-ranged. This situation characterizes order parameters of a thermodynam- 
ical system at the critical point of a second order phase transition (e.g., the gas-liquid 
transition at the critical temperature and pressure). For this reason, we call this class 
critical systems. 



Finally, if 

yd^xC(x) = s(0) = 0, (9) 

it is possible to show that 

a^{R)r^R'^, (10) 

with a < d. In particular, it is possible to show that in any case d — 1 < a < i.e. 
cr^(-R) ~ R''-^^ is the minimal scaling behavior for the number fluctuations versus R 
for any point process (all these considerations can be directly extended to include also 
any "genuine" continuous stochastic mass density field In this case, there is an 
exact balance between positive and negative correlations in the density fluctuations in 
such a way to have Eq. Q. Therefore, infinite wavelength density fluctuations vanish, 
which imparts a degree of "order" even to stochastic point processes that satisfy ©. 
At sufficiently small we have 

s(k)~F, (11) 

with 7 > 0. It is possible to show that a and 7 are related in the following way: (i) 
If < 7 < 1, we have a = d — 7; (ii) If 7 > 1, then a = d—1 (the "proper" condition 
for superhomogeneity) . For 7 = 1 there will be logarithmic corrections. 

Since for the class of systems that satisfy the number fluctuations increase with the 
spatial scale slower than in a large class of correlated and uncorrelated point process 
(e.g., Poisson distribution), we call them super- homogeneous or hyperuniform point 
processes. Note that super-homogeneous point processes are at a type of "critical" 
point, but one in which the direct two-point correlation function 4] rather than the 
covariance C(x) is long-ranged. 



IV. THE ONE-DIMENSIONAL MODEL 



In order to construct a super-homogeneous point process suitable for a complete study, 
we begin with a one-dimensional regular lattice of points, i.e., a chain of point-particles with 
constant spatial separation (lattice constant) a (see Fig^. The microscopic density for such 
a regular point process is given by 

+00 

n{x)= H^-Ja), (12) 

j=-oo 



7 



FIG. 1: Schematic representation of the one-dimensional lattice with lattice constant a. 



where a > is the lattice spacing. Clearly, such a set is not spatially stationary, but 
only possesses discrete translational invariance. However, it is the one-dimensional super- 
homogeneous point process with the lowest number variance as a function of i? 0. The 
global average density of the system is simply uq = 1/a. 

In order to obtain a stochastic super-homogeneous one-dimensional point process suitable 
for our study, we shuffle the lattice by introducing a random displacement field. That is, 
we move each point from its initial lattice position through a random displacement with a 
given probability distribution, each point being displaced independently of the others. In 
practice, if the initial position of the m^^ point is ma, the final one will be Xm = rna + Um, 
where random variable extracted from the probability density function (PDF) p{u). 

Note that the average density Uq is not changed by the application of the displacements, 
since the number of points in the system is conserved. 

It is possible to show (see Appendix |X} that if each point of a general initial spatial 
distribution is displaced from its initial position independently of the others with a PDF 
p{u), then the new power spectrum s{k) will be given by 

sik) = no{i-\p{m+simpm', (13) 

where si{k) is the initial power spectrum before the displacements and 

P{k) = / dup{u)e-'^'' (14) 



is the Fourier transform of p{u), i.e. the so-called characteristic function of the random- 
displacement PDF (for a more general discussion of the effect of a stochastic displacement 
field with arbitrary spatial correlation on a given point process see ^|). In general, we take 
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p{u) to be symmetric, i.e., p{u) = p{—u). Note that for all possible p{u), we have the limit 
condition p(0) = 1 and that for small k in the symmetric case 

p{k) ~ 1 - (15) 



with a = 2 and A = ^ if is finite, and where f{u) = dup{u)f{u) means the average 
over the uncorrelated displacements. For = 1, this is the case if p{u) decreases faster than 
for large Otherwise, if = +oo, i.e., p{u) ^ B\u\~^ for large |m| with 1 < (3 < 3, 
then a = (3 — 1 and jisf 

r+oo 

A = 2B rfxx-^(l -cosx), (16) 





where 5 is a positive constant. 

In the case of a lattice, it is well known |3| and simple to show that 

where the sum is over all of the integers m, except m = 0. Therefore, from Eq. |T^. the 
power spectrum of the "shuffled" lattice is 

_ /2'7rm\ 



a V a J \ a J 



(17) 



Recall that superhomogeneity (or hyperuniformity) of the point process is given by only the 
behavior of s{k) in the vicinity of /c = 0. Therefore, since in the first Brillouin zone the 
power spectrum of a lattice is identically zero (i.e. the first Bragg peaks are at \k\ = 27r/a), 
the small k behavior of s{k) is determined only by that of p{k). In particular, for \k\ < 27r/a 
and \k\ <^ (1/^) =^/° [cf. Eqs. fjl5|) and (fT^ J we have from the discussion above that 

2Ak'' {a = 2andA = f, if ^ < +oo 

s{k) = with <^ 2 _ (18) 

I a = /? — 1 and A from Eq. (fT6|) . if = +oo 

which always satisfies the superhomogeneity condition a > 0. In particular, for (3 > 2 we 
have 1 < a < 2, and the condition of minimal mass fluctuations-length scaling for point 
process in d dimensions (i.e. (t^{R) ~ R^'^) is satisfled (for (3 = 2 there are logarithmic 
corrections in L). 

In the case in which each point is completely randomly displaced inside its own unit cell, 
i.e., 

e(i-\u\ 
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i-1 



V . 



• — Q O • 



Q • 



^ X 

l+l 



i-3 



u. 



i-2 



U. 



i-1 



U. 



a 



FIG. 2: "Shuffled" lattice with the PDF of the uncorrelated displacements p{u) = ^ ^ — - with 
A < a. The filled circles represent the initial lattice configuration (i.e., a lattice with a lattice 
constant a), while the empty circles are the new positions of the points after the displacements Ui. 
The quantity Vi is the size of the final Voronoi cell of the point initially at the lattice position a ■ i. 

where 6{x) is the usual Heaviside step function, the final point distribution is not only super- 
homogeneous, but also completely statistically stationary (i.e., with a complete statistical 
translational invariance), even though the original lattice array was not. 

We first analyze the behavior of the fluctuations associated with the volumes of the 
Voronoi cells in the simple case in which (see Fig|21) 



The statistics of the Voronoi cells are relatively simple because no point is allowed to 
move into the unit cell centered at the initial position of another point. In what follows, 
starting from the results for this model, we will extend some of the results to the most 
general class of super-homogeneous point processes in any dimension. 

V. VORONOI-CELL STATISTICS 

As stated above, we start from a regular lattice of points with microscopic density given by 
Eq. ()12|). and displace each point independently of the others by applying to it a displacement 
whose PDF p{u) is given by Eq. (jl9|) . Taking the Fourier transform of this PDF yields the 
characteristic function p{k) to be exactly given by 




(19) 





fcA 
2 
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Consequently, applying Eq. (fTTj) we obtain 

/ 2 fkA' 



s{k)= - 
a 



—— sm — 
kA \ 2 



+ (20) 



27r ^ ^ /, 27rm\ / a /rrnrA 



j„ V a J VrnvrA 

m^O 



sm 



We can verify directly that, since sin(m7r) = for any integer m, only if A = a the contribu- 
tion to Eq. ()2()|1 coming from the Bragg peaks of the underlying lattice structure completely 
vanishes. In fact, it is the only case in which the point process is fully translationally 
invariant. 

We can now proceed to the evaluation of the statistics of the Voronoi cells. For a point 
process in any dimension, the Voronoi cell associated with a given point consists of the 
region of space closer to this point than to any other point. The collection of all of the 
Voronoi cells that tiles the space is referred to as a Voronoi tessellation. Clearly, in the 
initial lattice configuration, the Voronoi cell associated with each point coincides with the 
unit cell of size (length) a around each point. According to Eq. (fT^ . a randomly displaced 
point that was at the original lattice position ja (integer j) remains within its original unit 
cell. Consequently, we will always refer to this as point j. The size of its new Voronoi cell 
Vj will be given, by definition, by the size of the line segment that joins the point that lies 
exactly midway between the points j + 1 and j and the point lies exactly midway between 
the points j and j — I, i.e., 

.. = a + ^^^±l^, (21) 

where uj is the displacement applied to point j. The PDF fi{v) characterizing the size of 
the single Voronoi cell is formally given by 

fi{v) = J j ^ dxdyp{x)p{y)6 iv-a —j . (22) 

Use of Eq. ^ yields (see FigO)) 

if It; - a| > ■ 



2 



2(v-a) + A if-^<v-a<0 (23) 



2 

-2(f -a) + AifO<t;-a<f 



Let (...) denote the average over the realizations of the displacement field. Since we start 
from a deterministic point distribution (i.e., a lattice), this average is equivalent to the 
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FIG. 3: Representation of the one-cell size PDF /i(f) for our model. 



ensemble average over the final point process. In general, when also the initial state is a 
realization of a stochastic point process, the ensemble average over the final configurations 



by the double average ((•••)) Kiust be taken, where (...) is the average over the realizations 



of the initial point process, and (...) is the average over the displacements conditioned to 
the initial configuration. If the realization of the displacement field, seen as a continuous 
stochastic field with a value m(x) in each spatial point, is independent of the realization of 
the initial point distribution, the order of the two averages is totally arbitrary. It is only 
under this hypothesis that Eq. (fT^ is valid. 

Clearly, the average size of a Voronoi cell is given by 



Note that, as only finite up to A/2 jumps are permitted, only finite fluctuations for v 
are possible. The interesting question of whether infinitely large cell-size fluctuations are 
permitted in a super-homogeneous point process will be tackled in the next section together 
with other important aspects of Voronoi cells fluctuations. 

In the rest of this section, we analyze the joint probability distribution of two different 
Voronoi cells. In particular, we find an important "conservation law" for cell-cell correlations. 

In order to find the two-cell joint PDF f2{vi,Vj), it is important to note that in light 
of Eq. (j^TJ) Vi and Vj are two dependent variables only if |i — j| = 2. This means that for 





oo 



dvv fi{v) = a . 



The variance of the size of the Voronoi cell is given by 
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|i — j'l 7^ 2, we have 

f2{Vi,Vj) = fl{Vi)fi{Vj) . 

For j = i + 2, the PDF /2(fi, will be given by the integral 



f2iVi,V, 



i+2) 



dui^idui+idui+^p{ui^i)p{ui+i)p{ui+-i) 



5\Vi-a ^ j (5 \Vi+2 -a ^ J , (24) 

where is still given by Eq. ()19p. By performing explicitly the calculations and calling 
Wj = Vj — a for all j, it is simple to show that 



f2{Vi,Vi+2) = 



A - 2{wi + Wi+2) in 
A - 2wi 
A - 2wi+2 
A + 2{wi + Wi+2) in 
A + 2wi 
A - 2wi+2 




m 
in 



m 
in 



A2 
As 
A, 
A, 
A. 



elsewhere 



where (see FigH)) the Ai are the joint conditions: 



Ai = {wi>0 and < Wi+2 < -Wi + A/2} 

A2 = {0<Wi< A/2 and - Wi < Wi+2 < 0} 

^3 = {wi > and — A/2 < Wi^2 < —Wi} 

Ai = {wi<0 and - Wi - A/2 < Wi+2 < 0} 

A5 = {-A/2 <Wi<0 and < Wi+2 < -wi} 

Aq = {wi < and — Wi < Wi+2 < ^/2}- 



(25) 



The most basic and important quantity characterizing correlations between the size of 
different Voronoi cells is given by the correlation matrix Cij defined by 

Cij = {vi - a){vj - a) , 

where in this case the average is taken by using /2(fj, Vj) [cf. Eq. (j^ ]. Clearly, Cu = v'^ — v'^. 
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FIG. 4: Regions of the plane{wi 



a,Wi+2 = Vi+2 - a) where the joint PDF /2(fi,'Vj+2) 7^ 0. 



By direct calculation we have 



Cij 



24 for 2 = j 

48 



(26) 



for z=j±2 
for^^j,j±2 

We see that different Voronoi cells are either anti-correlated or uncorrelated in such a way 
that 

E C., = , (27) 

j = -OD 

i.e., positive and negative correlations must balance so that the sum of Cij over j is exactly 
zero. Because of the strong resemblance with the basic property Eq. Q of all the super- 
homogeneous point processes in arbitrary d dimensions, we expect that Eq. (j27|) is a general 
property of all super-homogeneous point processes in any dimension. 

To show that this expectation is indeed true, consider a spatially stationary super- 
homogeneous point process in d dimensions with average density of points no > 0. For 
such a point process, we know that the variance in the number of points N[R) in a sphere 
of radius R for sufficiently large R satisfies the relation 

(28) 



N\R)) - {N{R)y ~ with d-l<a<d 
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We will focus our attention on a given sufficiently large subset S of volume V (e.g. a sphere 
or an ellipsoid) and consider the number of points contained within it. The average value of 
this number is {N{S)) = rioV. Let us call Vi the volume of the Voronoi cell associated with 
point i. Since the set of point-particles is countable, we can arbitrarily label and enumerate 
them. By definition, {vi) = l/riQ. 

Let us now study the fluctuations of the quantity 

7V(5) 

u{s) = E 

i=l 

under the condition of superhomogeneity. Its precise value for a single realization will 
fluctuate from its average value given by 

{U{S)) = {N{S)){v,) = V. (29) 

In light of Eq. , we can write 

(|f/(5)-K|2)~\/t, (30) 
where it is to notice that {a/d) < 1. But from Eq. ()29|1 . we can rewrite 

{\U{S)-V\') = I Y: w^wA^V"^, (31) 

where, as in the one-dimensional case, Wi = Vi — I/iiq. This equation with a < d (condition 
of superhomogeneity), together with the fact that N{S) grows proportionally to V and the 
supposed spatially stationarity of the stochastic point process, implies directly that in the 
limit of an infinite subset S ,we have 

,N{S) y 

yl^(J2^^^^)=J2C^3=(^, (32) 

where Cij = {wiWj) and the last sum is extended over all of the point j of the system in 
the infinite volume limit. This result can be shown rigorously by various techniques, but it 
is sufficiently self evident to avoid having to present the mathematical details. This result 
is valid for any Voronoi cell i. In fact, while the matrix = (wiWj) depends on the way 
we have enumerated the points, in the case of spatially stationary point process, the sum 
J2j Cij does not depend on the enumeration. This is a quite interesting aspect of relation 

(El. 
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Therefore, in addition to Eq. Q, we have found another "sum rule" that apphes to all spa- 
tially stationary super-homogeneous point processes. To check that non-super-homogeneous 
point processes do not generally satisfy Eq. (|32|) is a very simple task. In fact from Eq. (|31|1 
it is simple to see that if a > d Eq. (j32|) can not hold. 



VI. LARGE CELL-SIZE FLUCTUATIONS IN SUPER-HOMOGENEOUS POINT 
PROCESSES AND VOID DISTRIBUTION 

In the previous sections, we analyzed the main properties of the one- and two-point statis- 
tics of Voronoi cells for super-homogeneous point processes. We found an important sum rule 
involving that the sum along any line or column of the Voronoi cells correlation vanishes for 
any super-homogeneous point process. In this section, we tackle two more important ques- 
tions about super-homogeneous point processes: (1) Can there be infinitely large Voronoi 
cells, or, equivalently, infinitely large voids, for super-homogeneous point processes? (2) Is 
it possible to find a functional expression for void size distribution linking the probability of 
having a void of a certain size to the correlation properties of the super-homogeneous point 
process? We will see that the answers to both questions are in the affirmative. 

The first question is motivated by the following facts: 



All the commonly known super- 



lomogeneous point processes (lattices, quasi-crystals 



{igI . the one-component plasma ^,0], (72— invariant processes [4] etc.) in the infinite- 
volume limit have only finite Voronoi cells and spherical voids; 

By taking Eq. (jH)), for a general stochastic mass distribution (continuous or point-like), 
it is possible to show [3] that, if s(k) ~ A;" at small k, then the wave modes which 
contribute essentially to create mass (i.e. number in point processes) fiuctuations on 
large spatial scales R satisfy: fc~l/i?ifn<l (and therefore cr^{R) ~ R'^~"'), and 
k kf) independent of i? if n > 1 (and therefore cr'^{R) ~ R'^~^ for all n > 1). In 
particular, fco marks the departure from the small k behavior of s(k) to its crossover 
to the large k behavior; in general "shot-noise" behavior for a point process, and a 
rapid cut-off to zero for a continuous mass distribution. Therefore, one might surmise 
that, at least in the case n > 1, voids much larger than the inverse of this cut-off wave 
mode fco are not permitted at all. This certainly is the case for the one- dimensional 
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model presented in the previous section in which Voronoi cells larger than twice the 
original unit cell (i.e., the inverse of the average density) are not permitted. 

However, more generally, we will see here that even in the case of Eq. ()28|) with a = d — 1, 
there are super-homogeneous point processes for which we can find spherical (or spherical- 
like) voids (and therefore Voronoi cells) that are arbitrarily large. Moreover, and importantly, 
for the case of shuffled lattices, we will derive mathematical relations between the probability 
of applying large displacements and the probability of having a void of the same size. This 
will permit to us to formulate an ansatz for the characterization of the whole class of super- 
homogeneous point processes in terms of the void size distribution. 

With this aim, we start again from the one- dimensional regular lattice of the previous 
section with lattice constant a = 1 and microscopic density given by Eq. |T2|) . We then 
again apply to it again an uncorrelated displacement field, but now we choose p{u) with 
an unlimited tail. As already shown in the Eq. ()18|) of the previous section, the final point 
process is always super- homogeneous satisfying the condition s(0) = for all possible p{u). 
With the aim of simplicity but no loss of generality in the final result , we restrict the analysis 
to the case in which p{—u) = p{u). 

Let us take the segment [0, 2R\ (i.e. the one-dimensional sphere of radius R) with R 3> 
a = 1, and ask for the probability W{R) that after the application of the displacement field 
no point is contained in it. Clearly, W{R) can be identified also with the probability that 
a randomly chosen void has a radius larger than R. Therefore, uj{R) = — '^^l^^ gives the 
PDF of the size (i.e. radius) of the voids. 

Given a point-particle initially at the lattice position m, it is simple to show that the 
probability Wm{R) that after the displacement u, it will outside of the segment [0,2i?], 
independently of m, is 




(33) 



where 




(34) 



J X 



Note that because p{u) is integrable over all the space 




lim (j){x) = and lim (j){x) = 1 



(35) 
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in any case. Since W{R) is the probability that all of the points in the system are outside 
of the segment [0, 2R] after the displacements, we can write 



+00 



W{R) = II [I- + 0(-m + 2R)] . (36) 

m=— 00 

In this equation, we can distinguish between two multiplicative contributions by writing 

W{R) = Wi{R)W2{R). 

The former contribution Wi{R) is given by the points initially outside the segment [0,2/2], 
and the latter W2{R) by those initially inside it. We show that the large- i? behavior of 
W{R) is determined essentially by this second contribution. 



1. Let us consider the first contribution: 

w,iR) = n [1 

m<0,m>2R 



-m) + 0(-m + 2i?)] 



(37) 



Because of the discrete translational invariance of the initial configuration, Wi (R) can 
be rewritten as 

"+00 "1 2 

Wi{R) = II {1 - ^{m) + ^{m + 2R)) (38) 



.m=l 



exp 



n (1 -0(m) +0(m + 2i?)) 

+00 

2 ^ In (1 - 0(m) + 0(m + 2R)) 



. m=l 

For any finite value of R, the convergence properties of the series 

+00 



^ In (1 - 0(m) + 0(m + 2R)) 



(39) 



m=l 



are given by the large-m behavior of In (1 — 0(m) + 0(m + 2R)). Because of Eq. ()35|) 
we can say that for sufficiently large m 

In (1 - 0(m) + 0(m + 2R)) ^ -0(m) + 0(m + 2R) . 

At this point we have to distinguish two sub-cases: 
(1) The PDF p{u) is such that 



/+00 
du \u\ p{u) < +00 . 
-00 
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In this case 



lim x(j)(x) = . 
2:^+00 ^ ' 



This imphes that 

+00 

(f){m) < +00 



m=l 



and so I]m=i l^^l^ ~ do. Therefore, to lowest order in 1/R, we can neglect 

(f){m + 2i?) with respect to (f){m) in Eq. (jHHjl and write 

iyi(i2)=Po>0, (40) 

where po = exp 2 1]m=i (1 ~ • Corrections to Eq. (jlUj) vanish for i? — >• +00. 

(2) If the PDF p{u) is such that 

r+00 

du \u\ p{u) = +00 , 
i.e. if p{u) = Bu~^^^ with < /3 < 1 for sufficiently large m, then 

lim xchix) = +00 . 

This implies that 

+00 

(j){m) = +00 , 

m=l 

being (f){m) ~ ^m~^ for large m. At any rate, the convergence of Eq. (j39|) . for any 
finite R, is still ensured by the following observation. In the limit m ^ 2R, we can 
write 

In (1 - 0(m) + 0(m + 2R)) ^ -0(m) + 0(m + 2R) ~ -2B Rm'^'^ , 

which guarantees the convergence of Eq. ()39|) . This implies that for large i?, iyi(i?) 
will have this main behavior 

Wi{R) ~ exp [-a{(5)R^-^] , (41) 

where a{(3) > 0. In particular, for /5 = 1, we expect that Wi{R) goes to zero for 
R — s> +00 as a power law. 
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2. Let us now analyze the second contribution to Eq. 



W2{R)= n [1 - + 0(-m + 2i?)] = n [1 - - 2i?) + 0H] , 

0<m<2_R 0<n<2/J 

(42) 

where in the last step we have adopted the change of variable n = 2R — m. Note that 
n — 2R < and that for R — > +oo with n fixed (j )(n — 2R) 1. Using the symmetry 
property p{—u) = p{u) of the PDF of the jumps |l8|, we can write 

(f){n -2R) = 1- (f){2R - n) . 



Therefore, Eq. ()42|1 can be rewritten as 
W2{R)= n [0(2i? - n) + 0(n)] = exp 

0<n<2_R 



J2 ln(0(2i?-n) + 0(n)) 

0<n<2_R 



(43) 



In order to evaluate J2o<n<2R^^ {4>{'^R — n) + (f){n)) let us approximate the sum by an 
integral as follows: 

r2R 



^ ln(0(2i?-n) +0(n)) ~ / dx In {(f){2R - x) + (j){x)) (44) 

0<n<2R "^^ 

rR 

= 2 dx ln((j)(2R-x) + (l)(x)) 
Jo 

Since (pix) is a decreasing function of x, we can introduce a further approximation by 
developing the In in Taylor series to the first order in 0(2i? — x)/0(x) : 



rR 

J2 ln(0(2/?-n) + 0(n)) ^ 2 / dx 

Q<n<2R ''^ 



In Six) + 



(t){2R-x) 
(p{x) 



(45) 



In general, the contribution given by the term (j){2R~x) / (j){x) can be neglected for large 
R with respect to the first one. We will use this approximation to study some simple 
but important cases: (A) a power-law tailed p(m), and (B) a generalized-exponential 
tailed p{u). 

(A) Let us consider the case in which p{u) = B u^^^^ with /? > for sufficiently large 
u. In this case for sufficiently large R one obtains 

fR 

/ dx\n(t){x) ^ -l3R\Tii2R) . 
Jo 

This implies in the same limit of large R 



\n{(f){2R-n) + (l){n))^-2/3R\n{2R). 

0<n<2R 



(46) 
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Corrections to this approximation can be neglected for large R as they are of the same 
or lower order than R. Finally, we can write 



W2{R) ~ e 



-2/3_Rln(2_R) 



-2(3 R 



We see that for any P > the quantity decreases faster than an exponential 

exp{—AR), and therefore this is the main contribution to the behavior of W{R) for 
large R, i.e., 

For /3 — > 0, the linear corrections in R to Eq. dominate giving W{R) ~ exp{—AR). 
This is well understood by considering that for /? — the final configuration of the 
point distribution will no longer be super- homogeneous, but Poissonian for which it is 
well known that the size of voids is exponentially distributed 



at sufficiently 



(B) Let us consider now the case in which p{u) ~ c|M|^exp 
large \u\ (in particular |m| > uq) with mq > 0, a > and any p. Again, we can use the 
approximation 

J2 ln(0(2i?-n) + 0(n)) ~ 2 / dx\n(j){x). 

0<n<2R "^O 

It can be shown by different techniques that for large x 3> 



a 



(X) ~ c^x^-'^^^ 



exp 



X 



0{uq/x)) . 



Therefore, for asymptotically large R (i.e. 2R ^ max[l,Mo]), the dominating behavior 
will be 

2mo 



^ ln(0(2i?-n) + 0(n))=. 



which implies 



W2{R) ~ exp 



2uo f R\ 



a + 1 \UqJ 

Since in this case Wi{R) approximately does not depend on R for large R, as in the 
previous case, W2{R) determines the behavior of W{R), i.e., 

a+V 



W{R) ~ W2{R) ~ exp 



2uo /RV 
a + 1 kuqJ 
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From the analysis of these two examples, we expect that, if p{u) is a PDF with an unlimited 
tail such that (as power laws and generalized exponentials) at large R 

lnp(7i?) =7'lnp(i?)(l + o(l)) (47) 

with 7,7' two positive related constants of order 1, then the following relation holds 

W{R) ~ exp{AR\np{R)) , (48) 

where A is a suitable constant depending on the average density of points uq and on the 
details of p{u). In fact, this result can be generalized to any other p{u) with unlimited tail 

The extension of this result to higher dimensions, in which again a regular lattice is 
perturbed by an uncorrelated displacement field characterized by a PDF with an unlimited 
tail, is straightforward when the PDF of the rf— dimensional displacement factorizes into a 
product of the PDF's of the single components Pd(u) = Ylf=ip{ui). In this case, by following 
the same procedure for the one-dimensional case, one can find that, given a cube of large size 
2R, the probability that it becomes completely void after the application of the displacement 
field is 

W{R) ~ exp{AR'^\np{R)) . (49) 

We expect that the above relation, with a suitable A, is also valid if instead of taking a 
cube a size 2R we take a sufficiently compact volume (e.g., a spheroid) linear size R. The 
mathematical treatment in the case of isotropic displacements Pd{u) = Pd{u) is more difficult, 
but we expect qualitatively the same result. We give only a rough sketch of this treatment. 
Let us take a sphere of very large radius R and, as above, factorize the probability W{R) 
that after the application of the displacements it becomes empty into the product of the 
probability Wi{R) that all particles initially out of the sphere stay out and the probability 
W2{R) that all the particles initially in the sphere go out of it because of the displacements. 
As in the previous case, we expect that W2{R) is the dominating factor for what concerns 
the large R decreasing behavior of W{R). This can be seen through the following arguments. 
In order to evaluate Wi{R) at sufficiently large R, we approximate the probability that a 
point, initially at a distance between r and r + Ar from the center of the sphere with r ^ R 
and Ar <C r, will stay out of the sphere after the displacement, as 

l-Mr)^i?^ 
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with Qd the complete spherical angle in d dimensions. Now the number of these particles 
in the initial lattice configuration is around rioQaf^'^^^Ar. Therefore, by taking the product 
over the spherical shells of thickness Ar for radii greater than R, we can write 



r)^R' 
d 



In the given limits Eq. (jHn|l can be reapproximated as 



Wi{R) ~exp 



no-fR" / drr^-'p,{r) 
d Jr 



(50) 



(51) 



In complete analogy with the one-dimensional case, it is simple now to see that, if Pd{u) 
decays faster than u~'^'^ at large u, then Wi{R) at asymptotically large L converges to a 
positive constant < po < 1- Instead if Pd{u) ~ Bu^'^^^ at large u with Q < (3 < d then 



Wi{R) ~ exp 



-aR 



where the constant a can be obtained approximatively by Eq. (j51|) . For what concerns the 
probability W2{R) we can say that for sure it must be smaller than the probability P{R) 
of the following event: all the particles (whose number N(R) is about ^{R/2Y) within a 
distance R/2 from the center of the sphere make a displacement u larger than R/2. This 
probability P{R) is (in the large hypothesis Eq. (jUj)) roughly given by 



P{R) 



+ 00 _ \^(^) 

duu'^ ^PdM ~ 

R/2 J 



exp 



CR''\npd{R) 



with C > appropriate and depending on d, and where we have considered the fact that, by 
definition, Pd{u) decays faster the at large u. On the other hand W{R) must be larger 
than the probability Q{R) that all the particles in the sphere make a jump of size larger 
than 2R. By similar reasoning one can find that 

Q{R) - exp [DR'^lnpdiR) 

with D another suitable constant depending on d. Since Pd{R) decrease to zero at large R, 
this shows at the same time that the decaying behavior of the factor W2{R) prevails on the 
one of Wi{R), and that again W{R) must have the form given by Eq. 

We recall now that in general in a Poisson point process in arbitrary dimension and with 
average density no, the probability Wp(y) that a given volume V is found empty of points 
is given by 



Wp{V) 



-noV 



(52) 
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For point processes that are essentially Poisson with primarily positive and short-range 
correlations, due to the only short-range clusterization of points, we expect a similar relation 
for sufficiently large voids, but with no replaced by an appropriate smaller constant 0, 2^ . 
On the other hand, in "critical" point processes, because of the strong clusterization of points 
at all scales due to large-scale positive correlations, we expect a larger probability of finding 
large voids than in the Poisson one. Therefore, for super-homogeneous point processes 
generated by displacing the points of a d— dimensional regular lattice in an uncorrelated 
manner, the probability that a compact volume of sufficiently large linear size R decays 
with R faster than in any non-super-homogeneous point process. 

This observation suggests the following general heuristic conclusion: 

• A point process is super-homogeneous if and only if its void size distribution W{R) 
satisfies the limit condition 

lim ^^ = 0. 

Moreover, the above discussion about the void distribution generated in a lattice by an 
uncorrelated but power-law-tailed displacement PDF suggests a second general heuristic 
conclusion: 



A point process for which 



\nWiR) 
lim ^ ' = 



not only is super-homogeneous but its power spectrum satisfies s(k) ~ k'^ at sufficiently 
small k with n > 2. However, this proposition cannot be inverted. In fact, n > 2 is 
obtained also in the case of power-law-tailed Pd{u), but with a finite variance. 

VII. CONCLUSIONS 

Super-Homogeneous stochastic point processes, and more generally super-homogeneous 
mass stochastic density fields, are a very important mathematical models of many sys- 
tems not only in material science and condensed matter physics, but also in diverse fields 
such as cosmology. For instance, slightly perturbed crystal lattices, quasi-crystals [3| , one- 
component plasmas 0, particular glassy systems, strictly jammed stochastic hard spheres 
configurations can all be seen as super-homogeneous point density fields. Cosmological 
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models predict a spectrum of the primordial mass density perturbations of the Universe 
typical of super-homogeneous systems 0,0, and super-homogeneous point processes (typ- 
ically perturbed lattices or glass-like particle distributions) are used as initial conditions in 
n— body simulations to study the mass collapse and the structure (e.g, galaxies, and clusters 
of galaxies) formation problems during the history of the Universe. 

Usually a point process is recognized to be super-homogeneous by studying the scaling 
behavior of its number fluctuation (t{R) with respect to the distance at asymptotically large 
spatial scales (see Eq. (fTUj)). or by analyzing the spatial integral of the density covariance 
C(x) (see Eq. (P)) or equivalently the power spectrum s(k) at small wave numbers (see 
Eq. (HH)). 

In this paper, we have characterized super-homogeneous systems by studying the sta- 
tistical properties of the Voronoi cells and of void size distribution. It is an important 
achievement because the knowledge of the statistical properties of Voronoi cells is an im- 
portant issue in many subjects of disordered materials. This task has been accomplished 
mainly with the detailed study of the so-called one-dimensional "shuffled lattice" i.e., a 
regular chain of particles whose particle are randomly displaced from their lattice positions 
with no correlations between the displacements. Inspired by the achievements obtained for 
these systems, we have generalized the main results to the whole class of super-homogeneous 
point processes in arbitrary spatial dimension. 

The main results that we have obtained can be summarized as follows: 

• For a particular subclass of one- dimensional "shuffled lattices", one and two Voronoi 
cell statistics have been solved exactly; 

• The correlation matrix Cij of the Voronoi cells of any super-homogeneous point process 
satisfies a sum rule J2j Cij = 0, which is independent of the way in which the single 
Voronoi cells have been labeled. This is a very important relation because it is a 
special property of only super-homogeneous point processes. Indeed, this sum rule is 
the Voronoi-cell equivalent of Eq. Q, which is the definition of a super-homogeneous 
point process in terms of the covariance function. 

• In contrast to the conventional picture of super-homogeneous systems, we have shown 
that arbitrarily large Voronoi cells or voids are permitted in the super-homogeneous 
class. This is true despite the fact that super-homogeneous point processes possess 

25 



the slowest number (mass) fluctuations-length scaling relation possible for any point 
process. 

• For the most general one-dimensional shuflled lattice, we have found the asymptotic 
form of the void size distribution and its dependence on the "shuffling" statistics. 

• This result for void statistics has been extended to higher dimensions, and suggests the 
introduction of two heuristic conditions to identify and classify any super-homogeneous 
point process in terms of the asymptotic behavior of the void size distribution. 

This last result together with the sum rule about the correlation matrix of the Voronoi cells 
are the two most signiflcant achievements of this study. The present analysis and results 
open the possibility for new studies on even more complex morphological characterizations 
of super-homogeneous point processes. 
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APPENDIX A 

In this appendix we give a brief derivation of Eq. (fT^ . For a more general analysis of 
the effect of a stochastic displacement field on the power spectrum of a given point process 
see Ref. H- 

Let us call n/(x) = J2iLi ^{^ ~ ^i) the initial microscopic density of a given point process, 
defined on the line segment [— L/2, //2], where —L/2 < Xi < L/2 is the position of the i*'' 
point-particle of the system before the application of the displacement field. Ultimately, we 
will take the limit L oo. Let us also suppose we know the global average density rig = 
\imL_^oo N/ L and the power spectrum si{k) of such a point process as defined respectively 
by Eqs. (j2I) and (jSJ. We now apply to each point i, independently of the others, a stochastic 
displacement Ui extracted from the probability density function p{u). The new microscopic 
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density will be 

H^) = ^S{x - Xi - Ui) . 

i 

By definition, the new power spectrum s{k) will be given by 



1 /^'^ \ 
s{k) = lim - ( Ve-*''(^»-^j+"«-"j) ) - 2'Knl6{k) , (Al) 
L^+oo L \ . . / 



where (...) stands for the average over all the possible realizations of the displacement field for 
a given realization of the miiia/ point process, and (...) stands for the ensemble average over 
all the possible realizations of the initial point process. In our hypothesis the displacement 
field and the point process are considered statistically independent, and hence the two 
averages commute and they can be taken in an arbitrary order. We will first take the average 
over the displacements by separating the diagonal contribution from the non-diagonal one 
in the double sum of Eq. (jAljl : 



1,N 1,N 

^Q-ikixi-xj+Ui-uj) = jv+ |p(fc)|^ ^'e"*''^''^"^-''^ , 

i,j hi 

where p{k) is defined by Eq. (fT^ and Y^'ij means the sum over alH = 1, and j = 1, A^ 
with i ^ j. Therefore, we can rewrite Eq. (jAlj) as 

|2 /l,Af ' 



s{k) = lim 



^(i-ipwiV^(i: 



g ik(^Xi Xj') 



2'nnl5{k) , (A2) 



where we have added and subtracted the term ^ |p(^)|^ = Ya=i in order to complete 
the double sum. Equation (jlHj) is recovered by noticing that the following relations hold: 

p(0) = 1 
iim — = Uq 

L— >+oo ij 

1 /^'^ \ 

lim - ( V e-'^^^-'^^ ) - 2T:nl6{k) = si{k) . 

L^+oo L \ . . / 

The extension to higher spatial dimensions is obvious. 
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